GUT-Scale Primordial Black Holes: Consequences and Constraints 



Richard Anantua^, Richard Easther^, and John T. Gibhn, Jr^'^ 
^Department of Physics, Yale University, New Haven CT 06520 and 
'^Department of Physics and Astronomy, Bates College, 44 Campus Ave, Lewiston, ME O424O 

(Dated: August 21, 2009) 

A population of very light primordial black holes which evaporate before nucleosynthesis begins is 
unconstrained unless the decaying black holes leave stable relics. We show that gravitons Hawking 
radiated from these black holes would source a substantial stochastic background of high frequency 
gravititational waves (10^^ Hz or more) in the present universe. These black holes may lead to 
a transient period of matter dominated expansion. In this case the primordial universe could be 
temporarily dominated by large clusters of "Hawking stars" and the resulting gravitational wave 
spectrum is independent of the initial number density of primordial black holes. 



Primordial black holes (PBH) produced immediately 
after the big bang IlL |2[ can decay via the emission of 
Hawking radiation [3, The initial PBH population 
is determined by the primordial perturbation spectrum. 
Bounds on the PBH population constrain inflation and 
other early universe scenarios, which generate this spec- 
trum 0, [y, 0, Q. These constraints follow from the 
lack of evidence for the present-day existence of PBH, 
and because decaying black holes disrupt nucleosynthe- 
sis, recombination and reionization At formation, 
a PBH must be smaller than the Hubble horizon, and 
the amount of material inside the Hubble horizon - and 
the maximal mass and lifetime of a PBH - increases 
as the universe expands. Very light PBH decay com- 
pletely before nucleosynthesis, and are consequently un- 
constrained. A PBH radiates any and all particles whose 
rest mass is substantially less than its current tempera- 
ture, including gravitons. If the PBH radiate massive, 
long-lived particles one obtains tight bounds on their ini- 
tial population [l^, 11, 12, 13, 14 1 but these limits are 
contingent upon assumptions about particle physics and 
quantum gravity, and other radiated particles will reach 
thermal equilibrium, erasing any memory of their origin. 
However gravitons emitted as the black hole decays can- 
not equilibriate and will always survive until the present 
day, producing a stochastic background of gravitational 
waves. Further, for some parameter choices the early 
universe has a transient matter- dominated phase, during 
which large clusters of PBH can form. In this case the 
resulting gravitational wave spectrum is independent of 
the initial fraction of black holes. 

The mass-fraction of PBH is denoted ^bh- Initially 
^BH = < /3 < 1. We assume that the remaining 
matter consists of radiation. PBH form if 5p/ p > 10^'^ 
at very short scales, and a number of inflationary mod- 
els have this property [l^ and references therein). 
Beyond this threshold (3 rapidly approaches unity. An 
appreciable gravitational wave background is generated 
even if (3 is very small, so given a model which predicts 
the existence of any small PBH, the signal discussed here 
is generic. 

For simplicity, we assume a PBH population whose 



mass is equal to the energy contained inside the Hubble 
volume at the instant they collapse. Recalling that = 
87rp/3Mp, and deflning p = Ef^^^, 



BH 



p 



327r . 



(1) 



which is the mass contained inside a sphere of radius 
1/H. Including grey body factors, Tgi, a Schwarzschild 
black hole emits (massless) particles with momentum k, 
reducing its total energy as 
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where h{s) counts the helicity/polarization states of a 
particle with spin s [l7, 18 1. The second line is the pure 
black body expression and g is the effective number of 
(bosonic) degrees of freedom, after grey body corrections. 
These suppress emission at larger s so g depends on both 
the mix of spin-states and total number of light degrees 
of freedom. For each state with s — 0, 1/2, 1, 2, the con- 
tribution to g is 7.18, 3.95, 1.62, 0.18 so graviton emission 
is an order of magnitude below a naive mode-counting es- 
timate. We ignore angular momentum, which enhances 
graviton emission [l9| . The physical wavenumber k is 
k/a{t), where k is the comoving wavenumber and a(t) is 
the scale factor. Integrating, 
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from which we can compute the lifetime 
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An upper bound on Ei^^a comes from the inflationary 
energy scale, which is constrained by the non-detection of 



2 




f [Hz] 



FIG. 1: Q.gn,{f)h^ with (from left to right) E-.^it = 
10^^ 10", 10^^ and 10^^ GeV, and /? = 0.001 and g = 1000. 
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FIG. 2: ng^{f)h^ with (from top to bottom) g = 10^ , 10^ 
lO'^ and 10^ In all cases /3 = 0.001 and Binit = 10^^ GeV. 
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FIG. 3: f7g™(/)/i^ with (from top to bottom) /3 = 10"^ 10"^ 
10"^, and 10"^^ The /3 = 10"^ and 10"*^ cases lie on top of 
each other. In all cases g = 10^ and -Bi„it = 10^^ GeV. 



Gravitational Wave Background: Denoting the number 
density of PBH by n{t), the energy density pbh = 
n{t)MBH{t)- We thus solve^ 
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along with equation ([2]). One obtains by an ap- 
propriate rescaUng. We finally compute the present -day 
spectral energy density of gravitational radiation 13, 23| . 
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where p is the overall density and pg^ is the energy den- 
sity in gravitational waves. ^ 

Figure [1] shows f2gto(/)ft.^ as a function of i?i„it. The 
gravitational wave power is substantial, and at very high 
frequencies. Roughly speaking, the temperature of the 
universe scales as l/a{t). A decaying black hole is much 
hotter than the surrounding universe, but the emitted 
gravitational waves are redshifted by the same factor as 
other radiation. Consequently, these gravitational waves 
have a higher frequency than the present-day CMB. Low- 
ering Ei„it increases the PBH lifetime, enhancing this dis- 
crepancy and pushing the gravitational wave signal to 
higher frequencies. The "dip" at very high frequencies 
arises because these quanta can only be sourced by a 
small black hole, and are produced in smaller numbers. 
Conversely, increasing g reduces the fraction of emission 
into gravitational waves, lowering Q,gw ■ As r is inversely 
proportional to g, the gravitational waves are emitted 
when the universe is smaller, increasing the subsequent 
redshift factor of the emitted radiation, lowering their 
present day frequency, as seen in Figure [21 

Early Matter Domination: The primordial universe is ra- 
diation dominated, whereas PBH scale like matter. Ini- 
tially, flBH oc a{t) until either Qbh ~ 1, or the PBH 



a primordial gravitational wave background in the CMB, 
which we (generously) take to be ~ lO^^GeV. At the 
lower end we are interested in black holes which decay 
prior to nucleosynthesis with time to spare for thermal- 
ization, so we need t < 100 s. With E,„,, = lO^^ Q^y, 
r w 29/(7 s, and the initial temperature is 18.8 TeV. Stan- 
dard model states alone give g ~ O(IO^) and in what 
follows we (conservatively) assume g > 10"^. Lowering 
Ei„it slightly ensures that the PBH will survive through 
nucleosynthesis, so we assume i?i„it > 10^^ GeV. 



^ Gravitons radiated by PBH are considered in [20l . I21II in a non- 
expanding universe. The former omits numerical factors, under- 
estimating the PBH lifetime and present day frequencies, and the 
latter focusses on present-day PBH. Also, [lf| considers gravita- 
tional waves generated by the formation of PBH. 

^ This quantity depends weakly on g*, the number of degrees of 
freedom after the universe rethermalizes. This differs from the 
g that fixes t, as a decaying PBH is much hotter than the sur- 
rounding universe. We take = 200, and plot Q,gw(f)h? , where 
h is the dimensionless Hubble parameter. 
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FIG. 4; The region in the {iJi„it,/3} for which a matter domi- 
nated period is aUowed is plotted for g = 1000 (white), along 
with the generic requirement that lO^'^ < -Ei„it < 10^". 




FIG. 5: Qbh and Qrad are plotted for a scenario with i5i„it 
10^^ GeV, P = 10"*, and g = 1000. 



reach the final phase of their evaporation and ^rad begins 
to grow. If the universe does become PBH dominated, 
all the radiation in the "late" universe will have been 
emitted by PBH, with the "original" radiation making 
a negligible contribution. In this case Qgu, = .86/(7 after 
evaporation. This is simply the fraction of the total emis- 
sion in gravitons and flg^^f) is independent of If /? is 
very small or g very large, the universe is always radia- 
tion dominated and fig^ < .36/(?. Setting prad ~ -E-itit; ^ 
matter dominated phase takes place if 



> 



15 Mf, 



(12) 



Recall that H"^ = l/At^ in a radiation dominated uni- 
verse. If the universe remains radiation dominated un- 
til the PBH have decayed, its grows by a(T + ii„it) « 
(^/tinit)^^^ with a(ti„it) = 1. Thus, rjgu,(/) decreases with 
P if the above inequality is not satisfied, as shown in Fig- 
ure m Figure H] shows the region of parameter space for 
which a matter dominated phase occurs, while Figure [5] 
shows VL bh and flrad for a specific scenario with a lengthy 



matter dominated phase. 

In a radiation dominated universe, H{t) ~ l/a(i)^. As 
always 1/H defines the physical Hubble scale while the 
comoving Hubble distance is a{t)/Hi„n. The number of 
PBH per initial Hubble volume is (3, so before matter 
domination, the number of PBH per Hubble volume is 
(3a{t)^. This number can be large: in FigureEJ /3 = 10~*, 
and a{t) — 10* before PBH domination, so there are 10^^ 
PBH within a single Hubble horizon. Perturbations grow 
in a matter dominated universe. A mode which is in- 
side the horizon and longer than the Jeans length has 
amplitude 5 oc if , -q = J dt/a{t) During matter 

domination, 5 o,{t), and short scales become nonlin- 
ear. Moreover, in order to ensure the formation of PBH, 
the initial amplitude of the perturbations is presumably 
substantially larger than the canonical 10~^ found at as- 
trophysical scales. A PBH dominated phase may thus 
be accompanied by the growth of nonlinear structure at 
sub-horizon scales, leading to the formation of large clus- 
ters of PBH. This situation is reminiscent of the present 
universe, with the decaying, clustered PBH playing the 
role of "Hawking stars" . 

The possibility that PBH cause a transient matter 
dominated phase has been discussed previously (e.g. (20j ) 
and a universe dominated by decaying PBH is in ther- 
mal equilibrium and thus a potential site for baryoge- 
nesis 2^ 25, [2^. Crucially, the formation of nonlinear 
over-densities could dramatically enhance the interaction 
rates between black holes [23,|2a|. If two PBH merge, the 
resulting black hole lives roughly eight times longer than 
the parent objects. If a typical PBH survives until shortly 
before the onset of nucleosynthesis, a small population of 
longer lived black holes is potentially troublesome. Since 
the lifetime of the PBH depends very strongly on the 
initial energy, we see from equation ([6]) that a factor of 
10 in T can be eliminated by increasing E^^it by a factor 
10^/^ « 1.5 but the lower bound on E^^^t will only change 
substantially if many PBH coalesce into single objects. 

To put a crude lower bound on the merger rate, recall 
that our horizon-mass PBH, have a Schwarzschild radius 
Vg = 2M/Mp which is equal to the initial Hubble length, 
\/H. Assume that PBH separated by an initial comov- 
ing distance of crs will merge, where c is a number of 
order unity. In a comoving region of radius erg , we ex- 
pect to find ~ c^P PBH, so volumes with N PBH will 
be (c'^/?)(^~^^ rarer than volumes with just one PBH. 
Thus, if we reach a matter dominated phase the fraction 
of the universe composed of PBH with mass NMBH{ti„it), 
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Unless P is close to unity this initial 



merger phase will not yield a long-lived population of 
PBH. However, correlations in the initial distribution of 
PBH [28] or the formation of large, large nonlinear clus- 
ters of PBH could substantially enhance the merger rate. 

Discussion: We show that light PBH which evaporate 
before nucleosynthesis lead to a high frequency gravita- 
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tional wave background. At present, this is of theoretical 
interest, given that this background is at frequencies far 
beyond the sensitivity region of LIGO, or proposed space- 
based interferometers such as LISA, which are the most 
sensitive gravitational wave experiments currently in de- 
velopment. However, the existence of plausible high fre- 
quency backgrounds motivates the development of novel 
detector technologies. The spectral density of this back- 
ground is substantial, and may exceed that obtainable 
from phase transitions or bubble collisions 22, 2^. Fur- 
ther, a light PBH population can lead to a temporary pe- 
riod of matter domination before the onset of nucleosyn- 
thesis during which clusters of PBH could form, leading 
to a cold phase during which the primordial universe is 
dominated by clusters of "Hawking stars" . 

Gravitational waves generated during preheating or 
parametric resonance at the end of inflation have recently 
received considerable attention 
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PBH thus provide a further mechanism by which inflation 
- if it sources perturbations which lead to the formation 
of PBH - may generate a high frequency gravitational 
wave background. Very simple models of inflation do not 
yield PBH and thus have /3 = 0, but current bounds on 
the running of the spectral index a = dus/dlnk are com- 
patible with PBH production Q. Further, these bounds 
are obtained by extrapolating the full inflaton potential 
from the region traversed as astrophysical perturbations 
are generated. This is not valid for models where inflation 
ends abruptly, and these scenarios can lead to substan- 
tial PBH production although the simplest models 
of this form often predict Ug > 1, which is in conflict 
with current data. Consequently, we simply treat /3 as a 
free parameter, although it would be computable in any 
well-specified inflationary scenario. However note that it 
need not be large - even if /? < 10~^° one may still have 
a lengthy matter dominated phase, provided E^^^ is at 
the lower end of the allowed range. 

The analysis here contains a number of simplifying as- 
sumptions. However, these do not affect our basic con- 
clusion, that a high frequency gravitational wave back- 
ground generated by Hawking radiation is the only sig- 
nature of a quickly decaying PBH population which cer- 
tainly survives until the present epoch. 
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